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The self-energy E(k, w), the fundamental function which describes the effects of many-body interactions 
on an electron in a solid, is usually difficult to obtain directly from experimental data. In this paper, we 
show that by making certain reasonable assumptions, the self-energy can be directly determined from angle 
resolved photoemission data. We demonstrate this method on data for the high temperature superconductor 
Bi2Sr2C'aCu20s+ x (Bi22f2) in the normal, superconducting, and pseudogap phases. 
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I. INTRODUCTION 

The propagation of an electron in a many-body sys- 
tem is described by the Greens function, G(k, oj) = 
1/ [u — £k — S(k, cj)], where is the bare energy of the 
electron and the self-energy S(k,w) encapsulates the ef- 
fects of many-body interactions. A detailed knowledge of 
S(k, oj) is of critical importance in elucidating the micro- 
scopic physics of the system. If its k-dependence is not 
important, one can obtain information about X from a 
probe like tunneling, which measures the density of states 
given by a k-sum of the imaginary part of G. This was 
exploited to get a very detailed microscopic understand^ 
ing of strong-coupling electron-phonon superconductorsEJ 
like lead. In general, though, if £ depends on k, then 
momentum averaged probes cannot be used to extract 
the self-energy. 

The only truly k-resolved probe is angle resolved pho- 
toemission (ARPES). Under the assumption that the 
"sudden" approximation applies (that is, one can ignore 
the interaction of the photohole with the outgoing photo- 
electron), for quasi-2D systems (since the component of 
the momentum perpendicular to the surface is not con- 
served in the photoemission process) , and assuming only 
a single initial state (one "band" ) , th^g the photocurrent 
can be written in the following forr 



J(k,w) = C k 



5k 



du'A{k',oj')f{oj')R{oj-oj') + B (1) 



where Ck is an intensity prefactor (proportional to the 
square of the dipole matrix element between initial and 
final states), A = (— l/VjimG is the single particle spec- 
tral function, / is the Fermi function, and R a Gaussian 
energy resolution function (photon monochromator and 
detector). The sum Ylsk is over a small window in k- 
space due to the finite angular aperture of the detector. 
B is the background, which contains extrinsic effects such 
as inelastic scattering of the photoelectrons (secondaries) . 

In this paper, we exploit Eq. (1) to determine the elec- 
tron self-energy, and illustrate this for ARPES data on 



the high temperature superconductor Bi2212. In Section 
II, we introduce the methodology which is necessary to 
extract the self-energy from the data. In Section III, we 
discuss the issue of background subtraction. In Section 
IV, various results are presented for Bi2212 in the nor- 
mal, superconducting, and pseudogap phases. Finally, 
some concluding remarks are offered in Section V. 



II. METHODOLOGY 

Let us assume we know A. Given that, we can easily 
obtain S. A Kramers-Kronig transform of A will give us 
the real part of G 



ReG{oj) = P 



+ °° d oj<^l 



(2) 



where P denotes the principal part of the integral. Know- 
ing now both ImG and ReG, then £ can be directly read 
off from the definition of G. 



IrriTi = 
ReT, = lu — e — 



ImG 



{ReG) 2 + {ImG) 2 

ReG 
{ReG) 2 + {ImG) 2 



(3) 



To obtain ReG using Eq. (2) we need to know A for all 
energies. From ARPES, though, we only know the prod- 
uct of A and /. (While unoccupied states can be studied 
by inverse photoemission, its resolution at present is too 
poor to be useful for our purposes). This is not a limi- 
tation if an occupied k-state is being analyzed and one 
can either ignore the unoccupied weight or use a simple 
extrapolation for it (except that only ReT. + e is deter- 
mined). On the other hand, one is usually interested in k 
vectors near the Fermi surface. Therefore a key assump- 
tion will have to be made. We can implement our proce- 
dure if we make the assumption of particle-hole symme- 
try, A(ek,u>) = A{— ek, — oj), within the small k-window 
centered at k^. Then, A is obtained by exploiting the 
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identity A(e k , u)f(u) + A(-e k , -w)f{-uj) = A(e k ,w), 
which holds even in the presence of the energy resolu- 
tion integration in Eq. (1). Note, this can only be in- 
voked at kp, and was used in our past, work to remove 
the Fermi function from ARPES dataaB, where it was 
denoted as the symmetrization procedure (note that the 
"symmetrized" data will correspond to the raw data for 
u> <^ —2.2kT). Although the particle-hole symmetry 
assumption is reasonable for small \u>\ where it can be 
tested in the normal state by seeing whether the "sym- 
metrized" spectrum has a maximum at the Fermi energy 
{Ep), it will almost certainly fail for sufficiently large 
u! > 0. Nevertheless, since we only expect to derive £ 
for uj < 0, then the unoccupied spectral weight will af- 
fect the result only in two ways. The first is through 
the sum rule J dujA(uj) — 1 which must be used to elimi- 
nate the intensity prefactor Ck in Eq. (1). From Eq. (3), 
we see that violation of the sum rule will simply rescale 
ImTi, but not i?eS due to the uj — e factor. Our normal- 
ization, though, is equivalent to assuming nk F =0-5, and 
thus does not involve "symmetrized" data. The second 
influence comes from the Kramers-Kronig transformation 
in Eq. (2), which is a bigger problem. Fortunately, the 
contribution from large uj' > 0, for which our assumption 
is least valid, is suppressed by l/(w' — uj). Further, for 
kp, ek=0 and thus ReY. is not plagued by an unknown 
constant. 

Some comments should be made about using real data. 
Data noise is amplified in the Kramers-Kronig transfor- 
mation in Eq. (2), and it is desirable to filter the data. 
We have found for our purposes that a wavelet trans- 
form works excellently in this regard, in that it provides 
smoothed data without any distortion of intrinsic spec- 
tral features such as the. quasiparticle peak. We employ 
a 'de-noising' algorithms which transforms the data into 
the 'wavelet domain' using class 6 complex Daubcchies 
waveletsLJ. Then, all wavelet components with absolute 
values below a certain threshold are set to zero and the 
data are transformed back into the signal domain. The 
threshold is set at a level which removes all (or most) of 
the noise from the data. The advantage of using a wavelet 
transform, over e.g. a Fourier filter, comes from the lo- 
calised nature of the wavelets in the signal and wavelet 
domain, i.e. the removal of noise from one portion of the 
data has no effect on intrinsic features elsewhere. 

Moreover, it is desirable to obtain a self-energy which 
is not artificially broadened in uj due to energy resolution. 
This is handled by deconvoluting the energy resolution 
out of the data using a maximum entropy methods based 
on the 'Cambridge Algorithm 'El which we have found to 
be quite stable. Here, the entropy of a distribution is 
defined to be S = — ^2 n Pn hip™, where p n is the inten- 
sity at point n. The algorithm locates the solution with 
maximum entropy, subject to its being consistent with 
the data when convoluted with the experimental resolu- 
tion. Consistency testing is done using the \ 2 statistic 
(C), C = J2 n ( D ™ -E n * R) 2 1 On, where D n are the data, 
F n * R is the solution (F n ) convoluted with the resolu- 



tion function (R), and is the variance of datum n. 
Since a completely flat solution has maximum entropy, 
the algorithm selects the smoothest ('deconvoluted') so- 
lution consistent with the original data and should only 
generate structures that are demanded by the data, i.e. 
those which are above the noise. To minimize the ef- 
fects of the resolution, we use a high resolution data 
set (<7=7.5meV, FWHM=18meV) in the low binding en- 
ergy range, and combine this with a lower resolution data 
set ((7=15meV, FWHM=35meV) to extend the spectrum 
out to higher binding energy (this takes advantage of the 
fact that sharp spectral structures only appear at low 
binding energies). The effects of broadening due to the 
finite momentum window can be minimized by looking at 
regions of the Brillouin zone where the dispersion is weak, 
which is the case considered here (our momentum win- 
dow has a radius 0.0457r/a). This will be less of an issue 
when considering data from new high resolution detec- 
tors currently becoming available, where the momentum 
window can be smaller in area by a factor of 25 or more. 

III. BACKGROUND SUBTRACTION 

We now illustrate our method by using data from the 
high temperature superconductor Bi2212. We choose this 
material because of its obvious interest to the condensed 
matter physics community, the novel electronic phases 
it exhibits as a function of doping and temperature, its 
lack of dispersion along the c-axis-which justifies the 2D 
approximation implicit in Eq. ( 1 )li3 , and our own strong 
familiarity with its spectra. On the downside, there is 
the background issue implicit in Eq. (1). Looking at the 
ARPES spectra over a large binding energy range, we 
see that the nea.T-Ep spectral features of interest to us 
ride on top of a large background. Not only is it too 
large to be ascribed entirely to the incoherent part of the 
spectral function A, the ratio of the spectral peak to the 
background changes with photon energy implying that 
most if not all of the background is extrinsic. This is sup- 
ported by the fact that the magnitude of the background 
is sensitive to the photon incident angle and polarization. 
Moreover, for k vectors where the spectral peak has dis- 
persed through Ep, this background is still present. It 
is flat in energy and extends all the way to Ep above 
T c , but is gapped in the superconducting state. Taking 
all of the above facts into account, the likely source of 
the background is scattering from other k vectors out- 
side the nominal momentum window, probably due to 
surface roughness and/or the incommensurate nature of 
the Bi2212 superstructure. 

There are a number of potential ways in which to sub- 
tract this background. An ideal way if one is along a 
symmetry axis (seldom the case) is to subtract data from 
perpendicular photon polarizations so as to recover that 
part of the signal which obeys the appropriate dipole 
selection rules. In practice, one is usually limited to sub- 
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tracting data from two perpendicular k vectors since the 
polarization is fixed. Moreover, the finite diameter of 
the momentum window, possible sample alignment er- 
rors, and the enhancement of noise due to subtracting 
two data sets, limit the effectiveness of this method. An- 
other possibility is to subtract data from an unoccupied 
k vector under the assumption that it is all background. 
The obvious problem here, besides the above mentioned 
amplification of data noise due to subtracting two data 
sets, is the strong variation of the dipole matrix elements 
in the Brillouin zoncO that can also act to modulate the 
intensity of the background from one k vector to the next. 

Because of this, we have instead explored models which 
capture the essense of the observed background, in par- 
ticular a step-edge (flat) background rand a "Shirley" 
backgrouncO The latter is of the forrrJlj 

POO 

I(u) = P(u) + c sh / duj'P{uo') (4) 

J uj 

where I is the total intensity and P that due to primary 
electrons (thus, one solves for P by simple matrix in- 
version). Although the step-edge background looks like 
the ARPES intensity seen for unoccupied k states (hence 
its motivation), it has the disadvantage of having three 
adjustable parameters (its height, and the position and 
width of its leading edge). Despite the fact that the 
Shirley background is designed to model secondary emis- 
sion, an unlikely source of the backgroundtSEJ, it is sim- 
ilar to the step-edge background, has the advantage of 
only one adjustable parameter-jaiid has been used exten- 
sively in previous treatmentsOLil 

To implement the background subtraction, the high en- 
ergy tail of the data is fit to a constant plus a Lorentzian, 
and then csh in Eq. (4) is varied such that this con- 
stant becomes zero. This results in a smaller background 
than simply forcing the intensity to be all background 
beyond some energy. This is done for data up to 0.5 eV 
where a minimum is seen in the spectrum, since beyond 
this, the spectrum rises and thus the "tail" becomes com- 
pletely buried under emission associated with the main 
valence band. For energies beyond 0.5 eV, we assume 
this Lorentzian tail when performing the integral (to in- 
finity) in Eq. (2) (real, not fitted, data is used below 0.5 
eV, of course). The purpose of this procedure is to avoid 
artificially forcing ImT, to zero at some cut-off (a power 
law tail is not used because the resulting integral would 
not be convergent). Once this background is subtracted, 
then the data are symmetrized (by adding the data at 
positive and negative energies), normalized (by invoking 
the sum rule), then Kramers-Kronig transformed, which 
is done analytically by assuming A to be linear between 
data points (the Lorentzian tail beyond 0.5 eV has an 
analytic transform, of course). A similar procedure is 
used for the step-edge background. The height of the 
step is determined by fitting the high energy data to a 
constant plus a Lorentzian. In the superconducting state, 
the position and width of the step's leading edge is deter- 
mined by fitting the low energy data to a Fermi function 



(whose "chemical potential" is the position and whose 
"temperature" is the width of the step) plus a Gaussian 
(modeling the spectral peak). In the normal state, the 
step-edge background simply reverts to a constant in the 
symmetrized data and so no low energy modeling is nec- 
essary. 

IV. RESULTS 
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FIG. 1. Symmetrized spectrum for overdoped Bi2212 
(T C =72K) at T=80K at the (tt, 0) - (tt, tt) Fermi crossing, 
with the line a fit to a constant plus Lorentzian. For visual 
purposes, it is shown for ui > 0, though we expect reliable 
information only for uj < 0. This applies to all the figures. 

In Fig. 1, we show symmetrized data at the (7r,0) — 
(tt, 7r) Fermi crossing for a T C =72K overdoped sample at 
T=80K, and thus in the normal state. Note that the 
spectral peak is centered at zero energy, consistent with 
being at &f with the zero of energy at Ep . The line is 
a fit to a Lorentzian plus a constant (flat background), 
and is an excellent representation of the data (with a 
HWHM of 55 meV). This Lorentzian spectral shape at k 
is sufficiently broad to make the quasiparticle ill-defined, 
but may seem unusual given, the supposedly expected 
marginal Fermi liquid formll3. We have always found 
Lorentzian fits to work well in the vicinity of the (tt, 0) 
point in the normal statetU Moreover, in Bi2201, where 
the normal state can be accessed over a large temper- 
ature range, we again find equally good Lorentzian fits 
even at low temperatures. The difference from optical 
conductivity dataMI, which do indicate a marginal Fermi 
liquid form, may be resolved by noting that the region 
near (tt, 0) makes little contribution to the in-plane trans- 
port due to the flat dispersion. In fact, near the direc- 
tion (0, 0) — (tt, tt), a case has been made that a marginal 
Fermi liquid lineshape can adequately describe the data 
if a backgro und , subtraction similar to what we employ 
here is donetSl!3. This points to the possibility of a vari- 
ation of the momentum dependence of E along the Fermi 
surface, which our method can in principle explore with 
the advent of new detectors with improved momentum 
resolution. Finally, we note that if we restrict away from 
small energies, a constant plus a power law fits the data as 
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well as a constant plus a Lorentzian. Typically, the (neg- 
ative) power, a, is such that \a\ < 1 (smaller for smaller 
doping) , which would be consistent with a non Fermi liq- 
uid line shapeE3£j. The advantage of the Lorentzian is 
that it goes through all the data, not just the higher en- 
ergy part, though this may be fortuitous if part of the 
"background" turns out to be intrinsic (the power law fit 
has the potential advantage of a smaller constant back- 
ground than the Lorentzian fit). A power law tail would 
also be divergent in Eq. (2), and thus would have to be 
cut-off (how to do this is not clear, since the tail is buried 
under the main valence band emission). This issue will 
hopefully be resolved in the future by doing a detailed 
analysis of the spectra as a function of photon energy, 
photon incident angle, and polarization to determine how 
much of the "background" is truly extrinsic. Once this is 
achieved, a closer representation of the true self-energy 
can be obtained. 
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FIG. 2. (a) Symmetrized spectrum (smoothed and Shirley 
subtracted) for overdoped Bi2212 (T C =87K) at T=14K at 
(-7T, 0) with (dotted line) and without (solid line) resolution 
deconvolution. The resulting ImT, and ReE are shown in 
(b) and (c). The dashed line in (c) determines the condition 
ReT, = u. 



In Fig. 2(a), we show T=14K symmetrized data for 
a T C =87K overdoped sample at the (n, 0) point (data of 
Ref. ||l|). We note the important differences in this super- 
conducting state spectrum, compared with the normal 
state spectrum in Fig. 1, due to the opening of the super- 
conducting gap, with the appearance of a sharp quasipar- 
ticle peak displaced from Ep by the superconducting gap, 
followed by a spectral dip, then by a "hump"E3 at higher 
binding energies (where_the normal and superconducting 
state spectra coincide&Ej) . This unusual dip-hump struc- 
ture is only seen near (tt, 0). The resulting £ is shown 
in Fig. 2(b) and (c). At high binding energies, one ob- 
tains a constant ImE as expected from the Lorentzian 
behavior above T c in Fig. 1. Note the very large value 
(~ 300 meV), much larger than that implied by Fig. 1 
(this is verified by the normal state spectrum, which has 
a much larger HWHM than the normal state spectrum 
of Fig. 1). That is, the magnitude of ImT strongly in- 
creases with reduced doping. Near the spectral dip, ImT 
has a small peakdollowed by a sharp drop, which we had 
earlier infercedtil from the spectral shape guided by fits 
to the dataa. This behavior is expected if the electrons 
are interacting with a spectral distribution gapped by 
2A in the superconducting state together with a sharp 
collective mode inside the 2A gap. The current results 
fully confirm the collective mode explanation proposed 
in Refs. |l] and ||. 

Despite this sharp drop below 70 meV, ImT remains 
quite large at low frequencies. That is, the quasipar- 
ticle peak is not resolution limited. It's flat behavior 



(w 6 



j 7 ) between 20 and 60 meV is consistent with 



the T e dependence of the quasiparticle peak width noted 
in Ref. ||. Then, below 20 meV, there is a narrow spike 
in ImT. This is the imaginary part of the BCS self- 
energy, A 2 /((x> + i0 + ), which kills the normal state pole 
at w=0. The resulting l/u> divergence of the real part 
ReT,, which creates new poles at ±A=32meV, is easily 
seen in Fig. 2(c). This is followed by a strong peak in 
ReT, near the spectral dip energy, which follows from 
the Kramers-Kronig transformation of the sharp drop in 
ImT. The strong peak in ReT explains why the low en- 
ergy peak in A is so narrow despite the large value of 
ImT. The halfwidth of the spectral peak is given by 
r = ImT/Z where Z = 1 — dReT/du; (the inverse of 
the quasiparticle residue). In the vicinity of the spectral 
peak, Z is large (~9), giving a T of ~14 meV. We note, 
though, that T is still quite sizeable, and thus the peak 
is not resolution limited. As the peak is dispersionless 
near (tt, 0)1211, this width is unlikely to be due to momen- 
tum resolution, which was verified by simulation. One 
could ask if it were due to an improper energy resolution 
deconvolution. This is highly unlikely, which was also 
checked by simulation. For instance, if one fits the zero 
energy spike in ImT to a constant plus a Lorentzian, the 
resulting Lorentzian is extremely narrow (with a HWHM 
of2meV). 

It is crucial to understand the extent to which our re- 
sults for T depend upon the choice of various background 
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functions. In Fig. 3, we compare JmS (as in Fig. 2(b)) 
for three different background choices: Shirley, step-edge, 
and no subtraction at all (for the last case, the spectrum 
is simply chopped off at 0.5 eV binding energy, and thus 
no Lorentzian tail). It is reassuring that all three results 
arc qualitatively similar (at higher binding energies, the 
unsubtracted case decays to zero because of the cut-off). 
There are some interesting quantitative differences of the 
step-edge background from the other two, in particular 
the step-like drop in ImT is more pronounced (resulting 
in a much more pronounced peak in ReT). This behav- 
ior is not very sensitive to the choice of the leading edge 
position and width of the step-edge background, and the 
result is quantitatively close to the theory of Ref. ^3|. In 
all cases, ImT is quite large at low energies, consistent 
with a quasiparticle peak which is not resolution limited. 
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FIG. 3. ImT as in Fig. 2(b) (with resolution deconvolu- 
tion), but for three different background subtractions: Shirley 
(solid line), step-edge (dotted line), and no subtraction with 
a cut-off at 0.5eV (dashed line). 

We have also looked at data from a T C =85K under- 
doped Bi2212 sample (data of Ref. ||). Below T c we 
find behavior quite similar to that of Fig. 2. Of more 
interest in this case is the so-called pseudogap pi 
where a gap is seen in the spectral function above T c \ 
In Fig. 4(a), we show T=95K symmetrized data at the 
(ir, 0) — (tt,it) Fermi crossing. One again sees (Fig. 4(b)) 
a peak in ImT at w=0, but it is broadened relative to 
that of the superconducting state, and the correspond- 
ing divergence of ReT (Fig. 4(c)) is smeared out. Such 
behavior would be consistent with replacing the BCS self- 
energy A 2 /(lu + i0 + ) by A 2 /(uj + iT a ). We have recently 
shown that such a-.self-energy gives a good description 
of low energy datau, and can be motivated by consider- 
ing the presence of pair fluctuations above T c . In fact, 
the T of Fig. 4 looks remarkably similar to the simple 
form proposed in Ref. |^, even over a large binding en- 
ergy range. Note from Fig. 4 that although the equation 
u) — ReT(uj) = is still satisfied at \uj\ ~ A, ImT/Z is 
so large that the spectral peak is strongly broadened in 
contrast to the sharp peak seen below T c . Actually, to a 
good approximation, the spectral function is essentially 
the inverse of ImT in the range \u>\ <^ 2A. We can also 



contrast this case with data taken above T*, the tem- 
perature at which the pseudogap "disappears". In that 
case, the spectrum is featureless, and the peak in ImT is 
strongly broadened. As the doping increases, this peak 
in ImT disappears. Further doping causes a depression 
in ImT to develop around uj = 0, indicating a crossover 
to more Fermi liquid like behavior. 
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FIG. 4. (a) Symmetrized spectrum (smoothed and Shirley 
subtracted) for underdoped Bi2212 (T C =85K) at T=95K 
(pseudogap phase) at the (7r,0) — {tt,tv) Fermi crossing with 
(dotted line) and without (solid line) resolution deconvolu- 
tion. The resulting ImT, and ReT, are shown in (b) and (c). 
The dashed line in (c) determines the condition ReT — ui. 



V. CONCLUSIONS 

In conclusion, we have proposed a method for deter- 
mining the self-energy S(k, oj) from ARPES data. Al- 
though several important assumptions have to be made 
(particle- hole symmetry, background subtraction), the 
method has the advantage that one can directly deter- 
mine T, rather than attempt to guess it by fitting the 
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Given the wealth of information one can ob- 
tain, we expect this procedure to be very useful in elu- 
cidating the microscopic physics of solids, particularly 
low dimensional strongly correlated systems where many 
controversies exist. Specifically, we find a non-trivial 
frequency dependence of E in the superconducting and 
pseudogap phases of the high temperature cuprate su- 
perconductors, which puts strong constraints on the mi- 
croscopic theory for these materials. 



M. R. Norman et al, Phys. Rev. Lett. 79, 3506 (1997). 
D. S. Dessau et al, Phys. Rev. Lett. 66, 2160 (1991). 
M. R. Norman and H. Ding, Phys. Rev. B 57, R11089 
(1998). 

H. Ding et al, Nature 382, 51 (1996). 

A. G. Loeser et al, Science 273, 325 (1996). 

R. Claessen et al, Phys. Rev. Lett. 69, 808 (1992). 



ACKNOWLEDGMENTS 



We thank Yuri Vilk for a key suggestion. This work 
was supported by the U. S. Dept. of Energy, Basic En- 
ergy Sciences, under contract W-31-109-ENG-38, the Na- 
tional Science Foundation DMR 9624048, and DMR 91- 
20000 through the Science and Technology Center for 
Superconductivity. 



Present Address: Dept. of Physics, Boston College, Chest- 
nut Hill, MA 02467 

J. R. Schrieffer, Theory of Superconductivity (W. A. Ben- 
jamin, New York, 1964). 

2 H. Ding et al, Phys. Rev. Lett. 74, 2784 (1995). 

3 M. Randeria et al, Phys. Rev. Lett. 74, 4951 (1995). 

4 M. R. Norman et al, Nature 392, 157 (1998). 

5 M. R. Norman, M. Randeria, H. Ding, and J. C. Cam- 
puzano, Phys. Rev. B 57, R11093 (1998). 

6 A. G. Major et al, J. Phys. Cond. Matter 9, 10293 (1997); 
Code supplied by A. G. Major. 

7 I. Daubechies, Commun. Pure Appl. Math. 41, 909 (1988). 

8 S. B. Dugdale, Ph.D. thesis, Univ. of Bristol, 1996; Code 
supplied by S. B. Dugdale. 

9 J. Skilling and S.F. Gull, in Maximum Entropy and 
Bayesian Methods in Inverse Problems, eds. C.R. Smith 
and W.T. Grandy Jr. (D.Reidel, New York, 1985), p. 83. 

10 H. Ding et al, Phys. Rev. Lett. 76, 1533 (1996). 

11 A. Bansil and M. Lindroos, J. Phys. Chem. Solids 59, 1879 
(1998). 

12 D. A. Shirley, Phys. Rev. B 5, 4709 (1972); S. Hufner, 
Photoelectron Spectroscopy (Springer- Verlag, Berlin, 1996). 

13 L. Z. Liu, R. O. Anderson, and J. W. Allen, J. Phys. Chem. 
Solids 52, 1473 (1991). 

14 M. R. Norman, M. Randeria, H. Ding, J. C. Campuzano, 
Phys. Rev. B 59, 11191 (1999). 

15 C. G. Olson et al, Phys. Rev. B 42, 381 (1990). 

16 C. M. Varma et al, Phys. Rev. Lett. 63, 1996 (1989). 

17 See also C. G. Olson et al, Science 245, 731 (1989). 

18 A. V. Puchkov, D. N. Basov, and T. Timusk, J. Phys. 
Cond. Matter 8, 10049 (1996). 

19 P. W. Anderson, The Theory of Superconductivity in the 
High T c Cuprates (Princeton Univ. Pr., Princeton, 1997). 

20 S. Misra et al, Phys. Rev. B 58, 8905 (1998). 



G 



